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The deviations S( m ("intermittency corrections") from classical ("K41") scal- 
ing ( m = m/3 of the m th moments of the velocity differences in high Reynolds 
number turbulence are calculated, extending a method to approximately solve the 
Navier-Stokes equation described earlier. We suggest to introduce the notion of 
scale resolved intermittency corrections S( m (p), because we find that these S( m (p) 
are large in the viscous subrange, moderate in the nonuniversal stirring subrange 
but, surprisingly, extremely small if not zero in the inertial subrange. If ISR inter- 
mittency corrections persisted in experiment up to the large Reynolds number limit, 
our calculation would show, that this could be due to the opening of phase space 
for larger wave vectors. In the higher order velocity moment (\u(p)\ m ) the crossover 
between inertial and viscous subrange is (10^m/2) _1 , thus the inertial subrange is 
smaller for higher moments. 
PACS: 47.25. Cg, 47.10. +g, 05.40. +j 



Experimentally turbulent flow has long been known to be intermittent 0. A 
signal is called intermittent, if there are relatively calm periods which are irregularly 
interrupted by strong turbulent bursts either in time or in space. Correspondingly, 
the probability density function (PDF) develops enhanced tails of large fluctuations 
and a center peak due to the abundance of calm periods, i.e., the PDF becomes 
of stretched exponential type instead of being Gaussian. This also means that the 
r-scaling exponents ( m of the velocity differences, 

< |v r (x)| m >=< |u(x + r,t)-u(x,t)| m >oc/ m , (1) 

do not vary linearly with m, namely as m/3, as was originally suggested by dimen- 
sional analysis of the universal, inertial subrange of fully developed turbulent flow 
0. Q]. Any deviations 5( m = ( m — m/3 are called intermittency corrections. 

Phenomenological intermittency models describe the measured intermittency cor- 
rections 5( m more or less successfully. For a detailed discussion see e.g. ||, j5|. But 
from our point of view an understanding of intermittency has to come from the 
Navier-Stokes equation. 

As full simulations for high Reynolds number (Re ~ 10 6 ) turbulence are out 
of range even for near future computers, one is thrown on approximations of the 
Navier-Stokes dynamics. 

The main idea of such an approximation has been introduced by us in [[5], |5|. 
Meanwhile we have considerably improved our approach and in this paper we employ 
it to determine the intermittency corrections 5( m . For completeness, we briefly 
repeat how our approximation scheme works. 

It starts from the common Fourier series in terms of plane waves exp (ip ■ x), 
p = (p^, pi = riiL -1 , rii = 0, ±1,±2, .... The periodicity volume is (27rL) 3 , L is 
the outer length scale. To deal feasibly with the many scales present in turbulent 
flow, we only admit a geometrically scaling subset K of wave vectors in the Fourier 
sum, K = U[Ki, thus u«(x, t) = J2 p ek u i(p, t) exp (ip ■ x). Therefore we have called 
our approximation scheme "Fourier- Weierstrass decomposition" [^|, |5|. K = {p^\ 
n = 1, A^} contains appropriately chosen wave vectors, which already have quite 
different lengths but dynamically interact to a good degree. The K\ = {p^ = 
n = 1, A^}, / = 1, l m axi are scaled replica of K which represent smaller 
and smaller eddies. l max is chosen large enough to guarantee that the amplitudes 
u(p l ™ ax , t) of the smallest eddies are practically zero. Of course, l max depends on 
the viscosity v and thus on Re. 

We solve the Navier-Stokes equation for incompressible flow (i.e., p ■ u(p) = 0) 
in the subspace defined by the wave vector set K, 

«*(P) = -iM ijk (p) M i (qi)M fc (q 2 ) - vp 2 Ui{p) + fi{p). (2) 

qi,q2 6lf,qi+q2=P 

The set Kq is chosen in a way that as many triadic Navier-Stokes interac- 
tions p = qi + q2 as possible are admitted. The degree of the nonlocality 
in p-space of any triadic interaction can be characterized by the quantity s := 
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max(p, qi, q 2 )/min(p, q 1 , q 2 ). 
nently we choose 



We allow for s up to 5.74. To force the flow perma- 
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e u(p,i) 



for p G Ki 
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as a deterministic, non-stochastic driving. Ki n C Kq only contains the wave vectors 
with the three smallest lengths. The corresponding amplitudes u(p, t) carry the 
largest energy. 

For the numerical calculations the times are measured in units of the largest 
eddies turnover time L 2 / 3 e -1 / 3 and the velocities in units of (Le) 1 / 3 . The Reynolds 
number can then be defined by Re = u^ 1 . The coupled set of (3 • (l max + 1) • N) 
equations (2) is integrated with the Burlirsch-Stoer integration scheme with adaptive 
stepsize. All averages are time averages, denoted by (...}. 

We remark that the density of the admitted wave vectors per p-interval decreases 
as 1/p in our reduced waveset approximation, whereas it increases as p 2 in full grid 
simulations, see Fig.l. But this shortcoming at the same time is the main advantage 
of our approximation, because many more scales than in full simulations can be 
taken into account. In || and we achieved Re = 2 • 10 6 , i.e., 3 decades. We used 
iV = 26, l max = 10 and N = 80, l max = 12, respectively. The main features of fully 
developed turbulence as chaotic signals, scaling, turbulent diffusion, etc., are well 
described within our approximation [|, HJ. In particular, our solutions show small 
scale intermittency. This is accounted for by a competition effect between turbulent 
energy transfer downscale and viscous dissipation ||. 

The main improvements of our new treatment in the present paper are: (i) The 
number of modes per Ki is considerably increased up to N=86 instead of N=38 in || 
or N=26 in || . Thus the number of contributing triades in eq. (2) is much larger, see 
table 1. (ii) We now also allow for nonlocal interactions in p-space. The nonlocality 
of the waveset K with respect to the Navier-Stokes dynamics (2) can be quantified 
by s max , defined as the maximum of the s- values of all contributing triades. In our 
former calculations |5j, |6| we had s max < 2, which means that eddies can at most 
decay in half size eddies, whereas now for s max « 6 (see table 1) sweeping of small 
eddies on larger ones (up to a factor of 6) is possible, (iii) We can now consider 
the individual u(p,t) instead of the whole shells = Y,p£K, u (p)> what had to be 
done in || || because of the smaller number of triadic interactions, (iv) We are now 
much beyond a shell model || || ||, since on the |p|-axis the elements of the wave 
vector subsets K\ interpenetrate and intermingle considerably. 

We now offer our results. 

The spectra (|u(p)| 2 ) and (|u(p)| 6 ) calculated with N=86 wave vectors in K 
and with Re = 125 000 are shown in Fig. 2a. For comparison the same spectra are 
also given for N=38 as used in our previous work (Fig.2b, 0). As expected the 
scatter becomes less with increasing N. We fit the spectra with the three parameter 
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functions 

(|u(p)| m ) = c m p~ u exp (-p/p D ,m)- (4) 

In Table 2 the fit parameters ( m and po,m are listed. The ansatz (4) is theoretically 
known || to hold for m=2. We find that it also holds for m > 2 with pD,m — ^Pd^I m 
as one can expect, if in the VSR the higher moments factorize. For the dissipative 
cutoff pd := po,2 we obtain p D = (lln)^ 1 , where n = (z/ 3 /e) 1//4 is the Kolmogorov 
length. This well agrees with the long known (experimentally [jlOj and theoretically 



TT| , |T2"|| ) crossover between the viscous subrange VSR and the inertial subrange ISR 
in the structure function D^ 2 \r) = (|u(x + r) — u(x)| 2 ) at r about IO77 . According 
to our finding pD, m = 2pD,2/ m the crossover pr>, m in higher order moments (\u(p) \ m ) 
occurs at smaller p, namely, approximately at (lOnm/2)^ 1 . The ISR for higher 
moments is thus definitely smaller. This does not necessarily mean that the ISR for 
higher order structure functions D m (r) is also smaller, because they are not simply 
connected with (\u(p)\ m ) via a Fourier transform as in the case m = 2. In passing by 
we remark that by properly renorming the wave vector p and the spectral intensity, 
the spectra (4) can be shown to the universal for all Reynolds numbers both in 
experiment JT3] and in full simulations and in our approximate Navier Stokes 
solution 0. 

The intermittency corrections 5£ m from our overall fit (4) are much smaller than 
the experimental ones around. At the other hand we observe here as in || that 
there is much intermittency in the signals, at least for small scales. We therefore 
determined the exponents ( m in (4) by fitting restricted p-ranges only. We suggest 
to introduce "local" ( m (p). These are defined by local fits of the type (4), using for 
each wave vector p the moments in the local p-decades \p/\/TO, p^/lO]. The cutoff 
wave vectors are kept fixed at their global values po = (H77) -1 , PD,m = 2pr,/m. 
Also, as before (see caption of table 2) we devide the local ( m {p) by Ca(p)-Q 

The astonishing results are shown in Fig. 3a: There are large intermittency cor- 
rections S( m (p) for the small scales (large p, VSR), only moderate intermittency 
corrections for the large scales (small p, stirring subrange SSR), but hardly any 
deviations for p in the ISR. 

The small scale intermittency is well understood [15, and was extensively 
discussed in ||. It is best seen in small scale quantities as for example in the 
energy dissipation rate e(x, t) or in the vorticity. Here we observe in addition that 
the intermittency corrections S( m (p) in the VSR remarkably well agree with the 
r-scaling exponents /i(m/3), defined by (e™/ 3 ) oc r~^ m ^\ which we had already 
calculated in |J. This means that for r in the VSR Kolmogorov's refined similarity 
hypothesis (RSH) v r oc (e r r) 1//3 seems to be well fulfilled. This also is in agreement 



with Kraichnan's flTTL Frisch's [18|, and our [19] objections against the RSH in the 



ISR, arguing that for r in the ISR v r is an ISR quantity, whereas e r still mainly is a 



1 Instead of denning S( m (p) = Cm{p)/C,z{p) — m/3 one could also take the deviation 5( m (p) of 
the Cm{p) from the linear behaviour as a measure of intermittency. It holds SCm(p) = Cm{p) — 
to C3(p)/3 = (3(p)3(m{p)- In the VSR it is Ca(p) ~ 1, so both definitions for the intermittency 
corrections essentially agree, but in the VSR we find Cs(p) > 1, so S( m (p) > 6( m (p). 
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VSR quantity. Thus a relation like the RSH should only be expected, if r is in the 
VSR and both v r and e r are VSR quantities. 



Our result is also consistent with latest full numerical simulations p0| , which 
find the RSH fulfilled. Note that in these simulations r is always in or at least near 
the VSR since Re is still small. And last not least our finding also agrees with the 



observation of Chen et.al. |2T[ that the RSH is less and less fulfilled the larger r 
becomes. For further comparison with experiment, see below. 

Before we interpret the behaviour of the S( m (p) in the stirring subrange SSR and 
in the ISR, we checked how our findings depend on various changes of our Navier- 
Stokes approximation: (i) To be sure that the SSR-intermittency does not depend on 
the kind of forcing (3), we compared with the alternative forcing f (p) oc u(p), again 
p G K in . We also took a random forcing, but the results did not change noticeably, 
(ii) We varied the set Ki n and allowed for more or for fewer modes which are stirred, 
but again there was no change, (iii) We varied the type of wave vectors in Kq 
and their number N as well as the maximal nonlocality s max of the contributing 
triadic interactions (see table 1). Again, no sizeable change. In particular, the 
intermittency corrections did not increase with increasing nonlocality of the triadic 
interactions as we speculated in |J. (iv) Different values of S( m (p) were only obtained 
when the flow field was not yet statistically stationary, see Fig. 4. In Fig. 4a we 
averaged over 7 large eddy turnover times only. The total rate of dissipated energy 
tdiss = ^Z)peA'P 2 (l u (p)| 2 ) still exceeds the constant input e by about 1%. In this 
case there are considerable intermittency corrections S( m (p) for all p, which go 
down drastically in the ISR if one averages over 70 large eddy turnover times, see 
Fig. 4b, where we had statistically stationary results. (Stationarity is identified from 
the balance between the total dissipation rate and the total input rate.) Similar 
observations have been made when analysing experimental signals [ 2^| . (v) To 
demonstrate how 8( m (p) varies from run to run we refer to Fig. 5. The deviations 
6( m (p) for p in the ISR are very small, but still seem to be significant, (vi) We 
decreased the degree of locality of the ( m (p) by fitting the larger range \pj V20, pV20] , 
see Fig. 3b. Again no qualitative change; S( m (p) now tends to become even smaller 
in the ISR. (vii) We artificially extended the ISR by putting v = and extracting 
the energy from the smallest eddies by using a phenomenological eddy viscosity 
as employed in ||. Now, as expected, 8( m (p) ~ also for the large p, i.e., the 
small scale intermittency really originates from the competition between transport 
downscale and the viscous damping, (viii) One might speculate that intermittency 
corrections in the ISR would occur if our Fourier- Weierstrass ansatz would not only 



be wave number but also space resolving as in [23]. But when doing this we found 



that the intermittency corrections observed in K3] vanish if the number N of wave 



vectors in K is increased Q. One should note that we include in fact some degree 
of position space localization since any Fourier representation with many modes 
already allows for localization in space. 

To have another check, we also calculated the scale dependent flatness F(p) = 
(|u(j9)| 4 )/(|u(p)| 2 ) 2 oc jr r C4+2C2 i jf there is intermittency, then 2( 2 > (4, thus F(p) 
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has to increase with p. In fact we find such an increase of F(p) in the SSR from 
F(p = 3) ~ 2.7 (< 3, a result achieved also in various full numerical simulations 



and experiments, see e.g. [14j for a recent reference.) to the value F(p) « 3.0 valid 
for a Gaussian distribution. For p in the ISR F(p) = 3 stays constant. Approaching 
the VSR by further increased p, the flatness now strongly grows 0. This can be 
understood as being due to the small scale intermittency, as we extensively reported 
in |J. This behaviour of F(p) well agrees with the above described findings for 

*C(p)- 

Finally we report how the flatness of the velocity derivative Fi t \ = 
((<9iWi) 4 ) / {(diUi) 2 ) behaves as a function of Re. For the large Reynolds numbers 
which we consider we find Fi t i = 3.15, independent of Re. This again means, we find 
no intermittency as models which are constructed to describe intermittency obtain 
an increase of Fx,i with the Re number in terms of the ISR intermittency exponent 
/i(2), F M oc Re^' 4 } see e.g. g§ 0. 

Two conclusions of our findings are possible. 

First, the very small if not missing intermittency might be due to our approxi- 
mation. Even in the present, considerably improved ansatz the larger wave vectors 
are still thinned out, cf.Fig.l. If this argument was valid, the ISR- intermittency 
would have been identified by us as an effect of the opening of the phase space for 
larger wave vectors. Consequently, there should be no intermittency in 2D turbu- 



lence, where the energy cascade is inverse - and infact Smith and Yakhot [25] do 
not find intermittency in numerical 2D turbulence. 

The second possible conclusion is that there indeed might be no intermittency 
in the pure ISR in the limit of large Re. Of course, if so that must be due to the 
particular form of the nonlinearity, namely the u • gra<iu-term in the Navier-Stokes 
equation. It provides energy transport both downscale and upscale which, as our 
solutions show, fluctuates wildly and with large amplitudes around a rather small 
mean value of downscale transport. This nearly symmetric down- and up- scale 
transport is broken on the large scales (i.e., in the SSR) due to the finite size of 
the system. The largest eddies do not get energy by turbulent transfer downscale 
but only deliver turbulent energy to smaller scales. The symmetry of transport is 
also broken for small scales by the competition with the viscous dissipation. May 
be that the symmetry breaking mechanisms cause the large and the small scale 
intermittency. Note that Galileian invariance is only broken by the boundaries, i.e., 
by the finite size of the system. Both dissipation and our forcing scheme keep it. 

This second possible conclusion is in agreement with a recent theory developed 



by Castaing et.al. ||26|| . This theory predicts that 5( m = in the limit of large 



Re, and that in this limit Fi ; i is independent of Re, which we find, as mentioned 
above. The value of the Fii-limit, if it exists, is probably larger than what we 
find. Vincent and Meneguzzi |27| calculate already F\^ = 5.9 for a Taylor Reynolds 
number Re\ ~ 100. Castaing et al. |26|] find from their data analysis, that the 
flatness F(p) increases as logF(p) oc (r/p) 13 with (3 oc 1/ log (Re\/75). Our finding 
F(p) m 3 in the ISR (for Re « 10 5 , Re\ w 9000) well agrees with the large Re 
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limit of this experimental behavior. Note that in experiment it is still (3 ~ 0.24 even 
for Re\ = 2720 [[^J. Also the measured intermittency correction S( m at least for 
m > 6 are not even for Re\ = 2720 |26|, ^8[. However it could well be that in 
experiment the intermittency corrections might be overestimated, because they will 
tend to increase if the averaging time is not large enough and the flow is not yet 



statistically stationary, see above, Fig. 4, and the remarks in Ref. |22[ . 

Our finding, that there might be three ranges for high Re turbulence - namely, 
the SSR with moderate intermittency, the ISR with practically no intermittency, and 
the VSR with strong intermittency - might also be supported by some experimental 
data arround. In Fig.6 the spectrum (|u(p)| 2 ) is shown, taken from Gagne's wind 
tunnel measurements [^9j with the very high Reynolds number Re\ = 2720. While 
in the ISR (2 = 0.67, i.e., 6(2 = 0, seems to be a good fit, for the large scales 
(SSR) the exponent (2 = 0.70 is more appropriate. In the VSR the exponential 
damping according to (4) is not separated, so that (2 cannot reliably be identified 
in that range. Note, that in the same experiment higher moments and the scale 
resolved flatness F(p), which is more sensitive to intermittency corrections, show 
intermittency also in the ISR. 



A similar interpretation seems possible by inspection of Praskovsky's |30[ data 
for (|v r | 6 ) measured at the also very high Re\ = 3200, see Fig.7: In the middle of 
the ISR we clearly have S( e (p) = 0, whereas in the VSR it takes the value 5( G = 0.31. 
This is precisely what one expects from the RSH, namely 5(q ~ \i ~ 0.30. In the 
SSR the intermittency correction is S(q(p) = 0.27, which is considerably larger than 
what we found in this range. May be this is due to the plumes, swirls, or other 
structures [31] which detach from the boundary in real flow and might increase the 
intermittency in the SSR. 

Finally, we remark that also the quasi-Lagrangian perturbation analysis of the 
Navier-Stokes equation, done by Belinicher and L'vov [|32|, leads to S( m (p) = for 
p in the ISR in the large Re limit. 

To summarize, it can not yet be ultimatively decided which of the discussed 
conclusions of our numerical data will turn out to be robust. Either there in fact 
is ISR-intermittency also for Re — > 00 as an effect of phase space opening for large 
wave vectors (which by construction of our approximation scheme we miss), or 
there is indeed no ISR-intermittency in the limit of large Re [26, [32| . To decide this 
alternative, it would be very helpful to at least allow some opening of phase space, 
e.g., to increase the number of wave vectors per level as log A; as already done in a 
2D approximate solution of the Navier Stokes equation [Q. If then intermittency 
does not show up again, we clearly have to favour the conclusion that there is no 
ISR-intermittency in the large Re limit as our results demonstrate. 
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Table 1 



N 


interaction 
number of 
triades 


^TYICLX 


b 


26 [HI 


39 


1.92 


300 


38 @, PJ 


102 


1.92 


170 


50 


273 


3.46 


80 


74 


741 


5.00 


70 


74 


729 


5.74 


70 


74 


858 


3.46 


70 


80 


783 


5.74 


70 


86 


966 


5.74 


65 



Characteristic parameters of several different wavesets K . N denotes the number 
of wave vectors in K Q . In the second column the number of triadic interactions 
p = qi + q 2 between the wave vectors of any one set K[ is given, I not too large 
or small to avoid edge effects. s max is the maximal nonlocality (definition see text) 
of the triadic interactions, b is the dimensionless constant in the structure function 
D{r) = (|u(x + r) — u(x)| 2 ) = 5(er) 2 / 3 , r G ISR. The experimental value is b=8.4 
[TO]. The larger values in our approximation can well be understood 0, the decrease 
of b with the increase of N is in keeping with that explanation. 



11 



Table 2 





(|u| 2 ) 

\ 1 1 / 


(|u| 3 ) 

\ 1 1 / 


(|u| 4 ) 

\ 1 1 / 


(|u| 6 ) 

\ 1 1 / 


(|u| 8 ) 

\ 1 1 / 


(lul 10 ) 

\ 1 1 / 




0.682 


1.021 


1.359 


2.034 


2.707 


3.370 


Cm/ C3 


0.668 


1.000 


1.331 


1.993 


2.651 


3.301 


5 Cm 


0.002 


0.000 


0.002 


-0.007 


-0.016 


-0.033 


5Cm(^62) 


0.022 


0.000 


-0.044 


-0.200 


-0.444 


-0.778 


5( m (ran - (3) 


0.038 


0.000 


-0.046 


-0.170 


-0.347 


-0.590 


PD,m 


597 


398 


299 


201 


152 


122 


2p D)2 /m 


597 


398 


299 


199 


149 


119 



Results from the fit (4) to the spectra (|u(p)| m ) obtained with N=86 wave vectors 
in K for moments up to m=10. v — 8 ■ 10~ 6 . The average is over 80 large eddy 
turnover times (skipping the first 100 turnovers). We fitted the p-range [0, 1000]. In 
general C3 is rather near, but not exactly 1 as it should be according to Kolmogorov's 
structure equation [0. We therefore calculate the intermittency corrections from 



renormalized exponents (m/Cs, namely, <5Cm = Cm/ (3 ~ wi/3. For comparison, the 
values for Kolmogorov's log normal model [|4j], <5Cm = — /xm(m — 3)/18, are also 
given, which are well known to fit the data for m < 10 with /i = /i 2 ~ 0.20. 
5( m (ran — (3) = — log 2 [l — x + x(l/2) 1_m / 3 ] are the intermittency corrections due 
to the random /3-model (x=0.125) [Q. pn,m is the dissipative cutoff, which agrees 
very well with 2pD i 2/m as shown in the last row. 
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Figure captions 

Figure 1 

Comparison of the subset of wave numbers admitted in the Fourier- Weierstrass 
decomposition (upper) with the complete p- spectrum (lower), either < |p| < 25 
or < |p| < 500. While the exact density of wave numbers increases as p 2 , the 
geometric scaling 2 1 Kq makes the density decrease, i.e., the smaller scales are less 
well resolved. 

Figure 2 

Spectra (|u(p)| m ) for m=2 (o) and m=6 (+). 

a) N=86 wave vectors in K , v = 8 • 10~ 6 , s max = 5.74, Re\ = 9030, averaging time 
from 100 to 180. The input set K in of the forcing (3) consists of the 12 shortest 
wave vectors. The lines are the fits (4). 

b) N=38, v = 8 • 10~ 6 , averaging time 150 large eddy turnovers, same forcing as in 
(a). 

Figure 3 

S( m (p) for m=2,4,6,8,10, bottom to top. Same data as in Fig.2a. The shaded ranges 
on the right show the Kolmogorov values [Q 5( m = —fim(m — 3)/18 for fi = 0.20 
through fi = 0.30, because in [^6j \i = 0.25 ±0.05 is given as "best estimate". In (a) 
the fit range is [p/V4~0, pvlO], in (b) the larger local range [p/y/20, py/20] is chosen. 

Figure 4 

5( m (p) for m=2,4,6,8,10 (bottom to top). N=80, v = 5 • 10~ 6 , s max = 5.74, Re x = 
13 550. The averaging times are (a) 7 and (b) 70 large eddy turnovers, respectively. 

Figure 5 

(a) S( e (p) and (b) S( 2 (p) for N=80, v — 5 ■ 10~ 6 , s max = 5.74 for different runs. For 
two runs the averaging time is about 30 large eddy turnovers, for another two runs 
it is about 70. The (weighed) means and the standard deviations are marked by a 
dia- mond and by error bars, respectively. 

Figure 6 

Spectrum E{r/rf) from [16] for Re\ = 2720. Three ranges can be identified. On 
the large scales the SSR with ( 2 + 1 = 1-70, i.e., some intermittency corrections, for 
moderate r the ISR with (2 + 1 = 5/3, i.e., no intermittency, and for small r the 
VSR. 

Figure 7 

Moments of the velocity differences (|iv| 6 )/r 2 (°) and of the energy dissipation rate 
( e r) (•) (arb. units) against r, taken from [27]. The slopes of these quantities are 
S(q(p) and fi = /12, respectively. Clearly, for (|tv| 6 )/r 2 there are three ranges VSR, 
ISR, and SSR, whereas the VSR quantity (e 2 .) does not show different ranges. The 
arrows correspond to 20r//L and 1/5, respectively. 
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